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The stochastic variational method is applied to particle systems and continuum 
mediums. As the brief review of this method, we first discuss the application to 
particle Lagrangians and derive a diffusion-type equation and the Schrodinger equa- 
tion with the minimum gauge coupling. We further extend the application of the 
stochastic variational method to Lagrangians of continuum mediums and show that 
the Navier-Stokes, Gross-Pitaevskii and generalized diffusion equations are derived. 
The correction term for the Navier-Stokes equation is also obtained in this method. 
We discuss the meaning of this correction by comparing with the diffusion equation. 

I. INTRODUCTION 

The variational principle plays fundamental roles in many branches in physics, such 
as analytic mechanics, field theory, statistical mechanics, quantum many-body problems, 
scattering theory, fluid dynamics and so on. Due to its elegant mathematical nature, the 
approach serves also as important guidance to formulate models, extracting the essential 
physical structure of the system in question. In particular, the method is indispensable 
when we deal with symmetries and the associated conservation laws of a system 1|. All the 
consequences of symmetries, at least in the classical level, are simply condensed in the fact 
that the Lagrangian (density) is invariant under the transformations of the corresponding 
symmetry. 

On the other hand, once some dissipative processes are involved, the usual variational 
approach looses its elegance. To include dissipative effects, one has to introduce some extra 
function, such as the so-called Rayleigh dissipation function or a time dependent external 
factor for a Lagrangian The choice of the functions are usually ambiguous. This modifi- 
cation is attributed to the fact that dissipation is essentially related to the energy conversion 
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from macroscopic degrees of freedom to internal microscopic motions, accompanying with 
the entropy production In this sense, dissipation is beyond the scope of the classical 
variational method jsl, 0] . 

Around third decades ago, a new variational method was proposed 5|,|6|], where dynamical 
variables are extended to stochastic quantities. In this stochastic variational method (SVM), 
the effects of the microscopic degrees of freedom to the macroscopic dynamical variables are 
represented in terms of noises, leading naturally to dissipation of the dynamical variables. 
The advantage of this method is that we do not need to modify Lagrangians by introducing 
extra functions, and, once the noise is specified for a given Lagrangian, the dissipative 
dynamics is derived in a straightforward way. 

Since then, SVM has been studied for several cases, exclusively for particle Lagrangians. 
It was shown that diffusion-type equations 7H9] and the Schrodinger equation 10Ml3| can 
be derived in the framework of SVM. In particular, in the latter case, SVM approach has 
been extensively discussed in relation to Schrodinger and Nelson's stochastic formulation of 
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quantum mechanics [14|, H5[ . 

On the other hand, it seems that, in spite of the interesting aspects of the method, the 
potential possibility of SVM has not been explored sufficiently Recently, we have shown 
;hat the method can be used to derive the Navier-Stokes equation for compressible fluids 



16l | . In the present work, we follow up the previous report, and develop the discussion 



of SVM in more general way. Starting from a short summary of the method for the case 
of particle Lagrangians, we extend the discussion to the case of continuum mediums. We 
show the detailed derivations of the Navier-Stokes, Gross-Pitaevskii and generalized diffusion 
equations in the framework of SVM. 

Moreover, SVM can calculate a correction term to the Navier-Stokes equation. The phys- 
ical meaning of the correction is discussed by comparing with the usual diffusion equation. 

This paper is organized as follows. In Sec. II, we summarize the classical variational 
method to be contrasted to SVM. In Sec. Ill, the application of SVM to particle systems 
of four different examples. One of the results obtained in Sec. Ill is equivalent to the 
Schrodinger equation. The same argument is applicable even when there is the gauge in- 
teraction as is discussed in Sec. IV. There the gauge fields are treat as the external fields. 
In Sec. V, SVM is applied to continuum mediums, and the Navier-Stokes, Gross-Pitaevskii 
and generalized diffusion equations are derived. The possible generalization of the Navier- 
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Stokes equation is further discussed by comparing with the diffusion equation. Section VI 
is devoted to concluding remarks. 

II. CLASSICAL VARIATIONAL METHOD 

To clarify the difference between the classical variational method and SVM, let us restate 
the classical variational method in the following form. Suppose that there exists a vector 
field f = f(r,t) where the velocity of a particle at the position x(t) is given by 

d 



dt 



x(f) = f(r,t)| 



r=x(t) 



(1) 



As is shown, this quantity coincides with the particle velocity at the position of the particle 
r = x(£), but is defined for all r at a given t. This is essential when we introduce the 
stochastic nature for the particle trajectory. Our strategy is to determine this field f by the 
variational principle. 

Let us consider the usual action for a particle, which is given by the subtraction of the 
kinetic term and the potential term, 



dt 



t a 



m ( dx(t) 
~2 I dt 



V(x(t)) 



(2) 



where m is the mass of a particle and V is the potential energy. 

Let us consider the variation of the position, x — > x + 5x, with the conditions 5x (t a ) 
<5x (tf,) = as usual. The variation of the action is then given by 
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dt 



ta 



(ix(t) / dx(i) 
m — - — o 



dt 



dt 



dt 



- m -f(x(t),t)-W(x(t)) 



5x, 



(3) 



and the stationarity condition of the action leads to 



m(d t + f(x(t), t) ■ V)f (x(t), t) + V V(x(t)) = 0. 



(4) 



Of course, if we introduce the following variable, 



v(t) = f(x(t),t), 



(5) 



i 



then Eq. (J4]), together with the definition of the particle velocity (PQ), can be casted into the 
usual form of Newton's equation of motion, 

dx(t) 



dt 



v(t), (6a) 



m Mt)_ + vr(x ( t )) = 0? (6b ) 

dt 

where we used the definition of the material derivative, 

^ = ^&f^ = W + f(x( t ), t ).V)f(x W , t ). (7) 

On the other hand, we can regard Eq.Q as a special value of the solution of a partial 
differential equation for the field f(r, £), 

m(d t + f (r, t) ■ V)f (r, t) + V7(r) = 0, (8) 

with some appropriate initial and boundary conditions. In the case of the classical variational 
principle, we do not need to calculate f(r, t) for all r at a given t to determine the classical 
particle trajectory. However, when we extend x(t) to a stochastic quantity, as we will see 
soon later, we need to know the values of f(r,t) not only for the trajectory of a particle, 
but also for general values of r. This situation remind us the path-integral formulation of 
quantum mechanics where the knowledge of action is required as for general values of r. 



III. STOCHASTIC VARIATIONAL METHOD FOR PARTICLE SYSTEMS 

In the classical variational method, the evolution of a particle is deterministic. In SVM, 
we consider the case where the trajectory of a particle is stochastic Q|. This is because 
we would like to approximately take into account the effects of some microscopic degrees 
of freedom, which do not appear as the dynamical variables in a given action, through the 
effect of noises. Therefore, instead of Eq. ([TJ), the time evolution of the particle is described 
by the following stochastic differential equation (SDE), 

dx(t) = f(x(t),t)dt + V2^dW t (dt > 0), (9) 

where the last term represents the noise. Properties of the noise will depend on the detailed 
nature of the underlying microscopic dynamics. However, if there are clear separation of 
microscopic and macroscopic time scales, we can assume that microscopic information is 
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approximately replaced by the Gaussian white noise as the most standard and simplest 
choice. Thus, dW t is differential of a Wiener process, satisfying 

E[dW t ] = 0, (10a) 
E[dW l t dWl] = 5 l3 dt. (10b) 

where E[ ] denotes the expectation value for the stochastic processes. The magnitude of 
the noise is characterized by v which is a constant. We call Eq.(j9]) the forward SDE, because 
it is defined only for dt > (see the discussion below). 

In stochastic systems, the particle distribution function is introduced as 

p(v,t) = E[5(v-x(t))}, (11) 

where the expectation value is taken over all the trajectories {x(t)} . Therefore, p{r,t) con- 
tains the information of the time evolution for a given initial distribution p(r, to)- Of course, 
the distribution of initial values {x(t )} for Eq. ([9]) should be given by p(r, t ) so that 
Eq.( |TT]) is valid at t = t . From this definition and using the properties of the noise, the 
evolution equation of p is expressed in terms of the following Fokker-Plank equation |l7l |. 

d t p(r, t) = V ■ (-f (r, t) + z/V)p(r, t). (12) 

As discussed in the previous section, the field f(r, t) is yet unknown and is determined 
by employing the variational principle. However, differently from the classical deterministic 
case, we now face to a new problem: the trajectories of stochastic variables are not smooth 
and are non-differentiable. Thus the time derivative is not uniquely defined. For example, 
we may consider the following two possible definitions for the velocity of the particle at t, 

v t = lim f ^, 13a 

v(f) = lim X(t + rft) ~ x(t) . (13b) 

Obviously when the particle trajectory x(t) is continuous and smooth, the two definitions 
coincide, v(t) = v(i), as is the case of the classical variational method. On the other hand, 
for the stochastic trajectory, v(£) 7^ v(i) in general. In spite of this, it was shown that the 
variational approach can be consistently formulated by completing the information of the 
backward SDE in addition to the forward SDE, as is shown in Ref. ||5] j^ . 
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Let us introduce a SDE for the time reversed evolution of x(£) as 



dx(t) = f(x(t),t)dt + V2udW t (dt < 0), 
where we assume the properties of the new noises dW t stay the same as before, 



E[dW t 



0. 



E[dWi dWl} = S ij \dt\. 



(14) 

(15a) 
(15b) 



We call Eq. (fT4|) the backward SDE. We can obtain the Fokker-Plank equation from this 
backward SDE similarly to the case of the forward SDE as 



d t p(T,t) = V.(-f(T,t)-vV)p(T,t). 



(16) 



If the above backward SDE, as a set of stochastic processes, describes the correct repre- 
sentation of the time-reversed process of the forward SDE, the Fokker-Plank equation f lT6|) 
should coincide with Eq. ( TT2|) . From this condition, we obtain 



f (r, t) = f (r, t) + 2iA7 lnp(r, t) + V x A(r, t) + B(t), 



(17) 



were A(r, t) and B(t) are arbitrary functions. In the absence of external forces and non- 
singular boundary conditions, we can simply omit A(r, t) and B(£), so that f (r, t) and f (r, t) 
are related as 

f(r,t) = f(r,t) + 2z/Vlnp(r,t). (18) 
In order to incorporate the forward and backward SDEs consistently in the variational 
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scheme of stochastic variables [5|, H5[, we introduce the mean forward derivative, 

"x(f + h)- x(t) 



LhcCt) = lim E 

/i-i-0+ 



and the mean backward derivative, 



Dx(t) = lim E 



h 



x(t) - x(t - h) 



(19) 



(20) 



instead of the two velocities, Eqs. ( jl3a|) and ( 113b|) . Here, E[F(t')\P t ] denotes the conditional 
average of a time sequence of stochastic variables {F (£') ,t a < t' < U} only for t' > t, fixing 
the values of F (£') for tf < t [15]. Thus, DF{t) is still a stochastic quantity, defined by the 
forward SDE. In fact, this definition leads to 



Dx(t) = f(x(t),t). 



(21) 
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Similarly, E[F(t')\J-'t] is the conditional average for the past sequence, defined by the back- 
ward SDE, fixing the future values, so that 

Dx(t) =f(x(t),t). (22) 

These velocities are used to replace the classical velocity in the kinetic term of the action. 
However, this replacement is not unique and each case corresponds to different physical 
scenarios. In the following, we discuss four different cases. 



A. CASE 1: Pure Forward Derivative 



Let us consider the case where the velocity dx./dt in the classical action is replaced only 
by the mean forward derivative Z)x(i). Then the action is given by 

"m 



dtE 



-Dx(t) ■ Dx(t) - V(x(t)) 



(23) 



Note that the products of the mean forward (backward) derivatives are independent of 
the choice of the discretization scheme such as the Ito, Stratonovich-Fisk and Hanggi- 
Klimontovich schemes by the definition of the mean forward (or backward) derivative 17|. 
Now we introduce the variation of the position of the particle as 



x(t) — > x A (t) = x(t) + Ag(x(t),t), 



(24) 



where A is an expansion parameter and g(r, t) is an arbitrary smooth function which is 
differentiable and satisfies the boundary condition, 



Substituting into the action (|23~]) ; we obtain 



(25) 



h = I + A / dtE [mDx{t) ■ Dg(x(t),t) - W(r)| r=xW • g(x(t),t)] + 0(A 2 



I + A 
I + X 



ta 



ta 



ta 



dtE 



dtE 



{-mD(Dx(t)) - W(r)| r=x(4) } • g(x(t), t)\ + 0(A 2 ) 
{-mDf(x(t),t) - W(r)| r=x(t) } • g(x(t),t)l + 0(A 2 ). 



(26) 



Here, we have used the stochastic partial integral formula la, ll8|] and Eq. (12~T1) . Using the 
Ito fomula (Ito, , mma) H, the fi tst te™ ,„ the in te,; can be calculated as 

Di (x(t), t) = (d t + f (r, t) ■ V - ^V 2 ) f (r, t) . (27) 



r=x(i) 



s 



Therefore, the stationary condition of the action under the variations leads to 

m[d t + f (r, t) ■ V - 2u(V lnp(r, t)) ■ V - uV 2 } f (r, t) = - W(r). (28) 

Here f(r,t) is replaced by f(r, t) using the consistency condition ( 118]) . Equation ( 128]) and 
the Fokker-Planck equation ( TT2]) constitute a system of partial differential equations for the 
unknown function f(r, t) and p(r,t). 

As stressed in the previous section, in the classical variational method, the final equation 
of motion requires only the velocity v(t) = f(x(t),t) and the knowledge of the field f(r,t) 
at arbitrary position r is not necessary In SVM, we, however, have to calculate the field 
f(r,t) itself. In fact, f(x(t),t) which appears in Eq. (128]) cannot be completely replaced 
by v(t) because of the second spatial derivative term, which did not appear in the classical 
variational method. That is, the fundamental variable is not v(t) but f(r,t) in SVM. The 
velocity on the trajectory is known only after solving the forward SDE with the solution 
ofEq. (EHJ). 

From Eq. (l2T]) . f (r, t) gives the mean forward velocity of the particle, obeying the forward 
SDE equation. Thus, one might expect that f(r, t) is the velocity field corresponding to the 
probability current. From the Fokker-Planck equation (Tl2l) . we identify such a probability 
current as 

J(r,t) = (f(r,t)-z/v>(r,t). (29) 

One can see that f(r, t) is not the velocity field parallel to the probability current, and 
neither f(r, t) is. However, the average of f(r, t) and f(r,t), 

f m (r,t) = f(r,t) ^ f(r,t) = f(r,t) - z/Vlnp(r,t) = f(r,t) + z/Vlnp(r,t), (30) 

is exactly related to the probability current as 

J(r,t)=p(r,t)f m (r,t). (31) 

That is, the velocity parallel to J(r,t) is given by f m (r,£). 
Equation (128]) is then re-expressed as 

p(d t + f m • V) /J - V-{vpdi f m ) - V-(r/pV /J) - V-[MV {y ^lnp)] = (32) 

m 

where is the i-th component of f m . In terms of the mean velocity field, f m , the Fokker- 
Planck equation for p is reduced simply to the continuity equation, 

% + V • (pf m ) = 0. (33) 
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Equations (1321) and (1331) constitute a closed set of differential equations for f m (r, t) and 
p(r,t). Note that Eg. ( 132]) contains the term which violates the time reversal symmetry and 
describes a dissipative process. 

B. CASE 2 : Pure Backward Derivative 

Another possibility is to use only the mean backward derivative to express the kinetic 
term, 



dtE 



771 ~ ~ 

-Dx(t)-Dx(t)-V(x(t)) 



Repeating the variational procedure analogous to the case 1, we obtain 



7 A = / + A/ dtE {-m£)f(x(t),t)-VV r (r)| r=xW }-g(x(t),t) + 0(Y 

Jta 1 J 

The stationarity condition of the action leads to, instead of Eq. fl28]) . 

m[dt + f (r, t) • V + 2i/(V lnp(r, *)) • V + z/V 2 ]f (r, t) = -W(r). 
Equation for f m becomes 



(34) 



(35) 



(36) 



p(d t + f m • V) Q + V-{updi f m ) + V-(z/pV / m l ) - V-[(i/p)V (z/ d t lnp)] 



?7i 



VV. (37) 



One can see that this equation describes the time reversed process of Eq. ( |32|) . In fact, Eq. 
( 1371) coincides with Eq. (1321) by exchanging the sign of ia 



CASE 3: Mixed Product of Two Derivatives 



a 



Hasegawa [2| studied the case where the kinetic term is written with both of Dx and Dx 



as, 



dtE 



m 



(38) 



Differently from the previous two cases, this kinetic term is symmetric for the exchange 
between D and D. The stochastic variation, then, leads to 



1 



(dt + f m • V)f l m + 2z/ 2 V i (p" 1/2 V 2 x /p) = W. 



m 



(39) 



Differently from Eqs. ( )32l) and ( l37j) . one can see that this equation is symmetric for the 
time reversed operation t -H- —t. This is the consequence of the symmetric replacement in 
the kinetic term. 
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One can confirm that Eq. (1351) is completely equivalent to the result of Ref. by noting 
the relation f m (r, i) = 2r/V5'(r, t). 



D. CASE 4: Average of Cases 1 and 2 - Schrodinger Equation 

As the last example, we consider that the kinetic term of the action is given by the 
average of the mean forward and backward derivatives [u|, 



dtE 



m £>x(t) ■ Dx(t) + Dx(t) • Dx(t) 



- V(x(t)) 



(40) 



In this case, the kinetic term is symmetric for the exchange of D and D, and hence the 
derived equation by the variation is symmetric for the time reversed operation. As a matter 
of fact, SVM leads to the following equation, 



Hi 



— \ dtf + (f • V) f + z/V 2 f + d t f + (f ■ V)f - uV 2 f \ = - VV. 



(41) 



This equation can be re-expressed as 



in 



(42) 



It is known that Eq. (142)) . together with the continuity equation ( 133|) . is (almost) equiva- 
lent to the Schrodinger equation [l0|. To show this, we introduce the following real function 
6, 



(43) 



2z/ 

Obviously, here we assumed that f m is irrotational, V • f m = 0. Then, from Eq. (l42p . we 
obtain, 



d t e + v (v^) 2 - ^(p^ 1 /^ 2 ^) + - — v 



2um 



0. 



This leads to 



d t e = -v iyef + ^(p-^v 2 ^) 



Ivra 



V + c(t), 



(44) 



(45) 



where c(t) is an arbitrary function of t, However, this is always absorbed into the definition of 
6. Thus without loss of generality, we can set c(t) = 0. On the other hand, the Fokker-Plank 
equation is reduced to 



d t In ^p = -2z/(V In y/p) ■ V6 - vV 2 9. 



(46) 
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Now we introduce the wave function as 



Then, from Eqs.( l45p and (146]) . the evolution equation of ip is given by, 

1 



(47) 



idti) — v 



2wm 



0. 



(48) 



By choosing v = H/(2m), we obtain the Schrodinger equation, 

2 



ihd t ip 



h 

2m V i 



V 



V 



in. 



(49) 



Note that, by the definition of i/j given by Eq. ( l4T|) . |^| 2 gives naturally the probability 
density of the particle. 

Of course, if the flow is rotational, we cannot introduce the phase of the wave function 
as Eq. (T43]) . Thus the equation obtained in SVM is not completely equivalent to the 
Schrodinger equation. 

The above derivation of the Schrodinger Equation in terms of SVM opens an interesting 
possibility for a new formulation of quantum mechanics, and extensive works have been 
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done so far |9|, HOJ , but in this present paper, we do not pursue further this direction. We 
only mention that there is an ambiguity for the choice of Lagrangian in the derivation of 
the Schrodinger equation. Let us consider the following action 

m ( £>x(t) • £>x(t) + £>x(i) • Dx(t) 
~2 \ 2 



dt 



a 



-(Dx(t) - Dx(t)) ■ (Dx(t) - Dx(t)) 



-V(x(f)) 



(50) 



where a is a constant. This action is reduced to the previous action ( 140]) in the vanishing 



a limit. As is shown in Ref. 



131 ] . the stochastic variation of this action again leads to the 



Schrodinger equation by choosing 



1 



h 



v 



(51) 



'I - a/2 2m 

Here, we discussed the wave function for the action given by Eq. ( 140 p . The wave function 
i[) defined by Eq. (I47p can be introduced even in the other cases discussed before. However, 
only this fourth case leads to the known form of the Schrodinger equation. For example, 
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when we apply the same argument to Eq. (|39|) . we obtain (again assuming that f m is 
irrotational) 



ihd t ip 



ft 2 



</>+ — 

2m 



1 1 



I 2 Iri 2|^l 4 



(V|V| 



2\2 



(52) 



Since this equation is not linear in ifi, the superposition principle required in quantum me- 
chanics is not satisfied. 

IV. PARTICLE INTERACTING WITH GAUGE FIELDS 



The gauge symmetry is an important principle to describe interactions among elementary 
particles. Thus whether SVM can be applied consistently keeping the gauge invariance is 
a crucial question. In the case of the Schrodinger equation, it is known that SVM leads to 
the correct wave equation with minimal coupling of a charged particle, starting from the 
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12|. 



corresponding classical action as sketched below 

Let us start with the action of a classical particle interacting with the external gauge 
fields, 



dt 



m 



,. _v . A - V - e0j , (53) 
where <p an d A are the scalar and vector potentials of the electromagnetic fields, and V is 



another potential energy, respectively. Here we used the CG 



check that this action is invariant for the gauge transform as 37|j 



Gauss unit. One can easily 



A — ► A' = A + V X , 

— > = 0- -^rX, 
c at 



(54a) 
(54b) 



where x is an arbitrary function. 

For the stochastic variations to derive the quantum mechanical wave equation, we replace 
this classical action by, 

m Dx • Dx + Dx ■ Z)x e £>x + Dx 



dt 



A-V 



(55) 



2 2 c 2 

The reason for this choice is that the action should be symmetric under the time reversal. 
Following the standard procedure, we find that the stochastic variation of the action becomes 

- / 1 + f 



I+X / dt 



2 v ; 2c 



-VA 1 - W - eV0 



?+0(A 2 ), 
(56) 
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so that we obtain 

(d t + f m ■ V)f m - 2 I /V(p~ 1/2 V 2 v /p) + -V(V + e<f>) = -—d t A + — f m x (V x A). (57) 

m cm cm 

The second term on the right hand side represents the Lorentz force. 

In order to derive the corresponding Schrodinger equation, we introduce the phase 9 by 

f m + ^A 



2u 



instead of Eq. ( 143]) . Then, the equation corresponding to Eq. (I45p becomes 



d t d 



2cmv 



A) +KP- 1/2 V 2 ^) 



2vm 



(V + e</>). 



On the other hand, the equation of In J~p is 



d t In J-p = -2z/(V In yfp) ■ (V9) - uV 2 6 + —(V In Jp) ■ A + — (V • A). 

cm 2cm 

Finally, we find the equation for the wave function if; defined in Eq. (I47p becomes 



ihd t ip 



1 

2m 



e \ 2 

-ihV - - A + V + et, 

c / 



(58) 



(59) 



(60) 



(61) 



where again, v = ft,/ (2m). This is exactly the Schrodinger equation of a charged particle, 
with the minimal coupling of the external gauge fields. 

Note that, in the derivation of the Schrodinger equation, we have assumed that the 
vector field f m is irrotational. On the other hand, in the derivation here, as is seen from the 
definition of the phase (|58|) . f m can be rotational but f m + eA/(cm) should be irrotational. 
This fact shows that the well-known Aharanov-Bohm effect is closely related to the rotational 
nature of the flow of probability density. 

V. CONTINUUM MEDIUM 



So far, we have discussed the application of SVM to particle Lagrangian. In this section, 
we employ to a system composed of a continuum medium [l6|. 

Let us consider a continuum medium described by the mass density p m (r, t) and the 
collective flow velocity field v(r, t). Then the classical action of this system (in the Euler 
coordinate) is given by a functional of the two fields, p m and v as 



t a 



dt / cf r 



Pm(r,t) 



v 2 (r,t) -e(p m (r,t)) 



(62) 
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where e is the potential energy (internal energy) density of the medium. For the sake of 
simplicity, we assume that the potential energy is a function only of p m . To perform the 
variational procedure directly from Eq. (162j) . we have to add the constraints for p m and v 
as is done in Refs. l|, |3| . Due to this, in order to introduce S VM, it is more convenient to 
use the Lagrangian coordinate system {R}, rather than the Euler coordinate system {r} to 
specify the elements of the continuum medium. As is well-known, the Lagrange coordinate 
system is attached to the medium in such a way that the coordinate R represents the initial 
position of a element. In terms of this Lagrangian coordinate system, the above action is 
expressed as a functional of the trajectory of every element, 

2 



/ [x(R;f)] = [dtj d 3 R Po (n) \ - J £(Pm(x(R;t)) 



(63) 



where x(R;t) represents the trajectory of the element specified by R, po(R) is the initial 
mass density and J is the Jacobian of the transformation, 

R ->■ r = x (R, t) , (64) 

that is, J = J(R;i) = det|<9r/<9R|. Assuming that there is no chaotic flows, the conservation 
of the mass of the system is expressed as 

d 3 r p m (r,t) = d 3 Rp (R), (65) 

so that we have the following relation, 

p m (r,t) = jPo(R). (66) 

It is easy to see that the classical variational principle for the action (163]) with respect to 



the trajectories x(R;i) leads to the Euler equation [16]. 

In terms of the Lagrangian coordinate, the implementation of SVM is straightforward. 
First the trajectory of the Lagrangian coordinate is determined by the following forward 
and backward SDEs, 

dx = f(x,t)dt + V2udW(t), (dt > 0), (67a) 
dx = f(x,f)df + V2udW(t), (dt < 0), (67b) 

where cfW(£) and cfW(t) are noise terms whose correlation properties are same as Eqs. ( TTUj) 
and (!I5]) . As before, we require that unknown fields f and f satisfy the consistency condition 
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(fl8"|) . Then the continuity equation obtained from these Fokker-Plank equations is expressed 
with the mass density as 



dp m 
dt 



V • {prrSm) = 0. 



(68) 



From now on, to simplify the notation, we will omit showing explicit the (R, t) dependence 
of the trajectories, x = x (R, t) , and the density, p m = p rn (R, t). 

Similarly to the previous section, we will replace the kinetic term with the forward and/or 
backward mean derivatives. For example, we chose 



la 



I,, 



As the variation, we consider 



dt / rf 3 Rp (R) 



1 



Dx. ■ Dx — J e(p Ti 



(69) 



x (R, t) — ► x A (R, t) = x (R, t) + Ag(x (R, t) , t), 



(70) 



where A is an expansion parameter and g(r, t) is an arbitrary function. Under this variation, 
the mass density p m is changed as, 

dg l 



S\p m = Pm (r A , t) - p m (r, t) = ^ A 



0(A 2 ). 



(71) 



Here Ay is the cofactor (ij) of the Jacobian J, which satisfies the following properties 

dJ 



d{dri/dR r 



(Vr,- x Vr fc ) r 







E dR~ Aij = * ' {frj x frk) = °' 



i * i 1 

j 3 



E E = E J<5 *^ = 



(72a) 
(72b) 

(72c) 
(72d) 



where V denotes the gradient with respect to R which should be distinguished from V. The 
indices k) in Eqs. (I72ap and (172bl) should be taken as cyclic. By using these properties, 
the action under the variation is expressed as 



/ - A / dt d 3 R p m0 E 



Ag ■ ( D£>x(R, t) + —V (-— 

V Pm V d ( l lPm)Pm 



0(A 2 ). 



(73) 
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Then the stationary condition for the variations 

DDx(R, t) + —V ) = 0, (74) 

p m V d{l/p m )p m J 

leads to (conf. [28]) 

Pm {d t + f(v,t)-V r -2v(V\np m (r,t))-V-vV 2 ] f(r,t) + v(- — j —) =0. (75) 
By using Eq. (I30j) . this can be re-expressed as 

Pm(d t + f m • V)P m - dj^pmdifj + v Pm dj / m l ) - ^ dj^iyp^djiydi lnp m )) 

= -V f-^-) . (76) 



d(l/p m ) p 

where p m is the solution of the continuity equation ( 168]) . as before. Because of the asymmetric 
entry of D and D in the action, the derived equation violates the time reversal symmetry. 



A. Navier-Stokes Equation 



The above scheme is applicable to derive the Navier-Stokes equation (NS). A derivation of 



NS for incompressible fluids in SVM has already been worked out in Refs. 



3,Q- 



Recently, 



in Ref . [16( , the present authors have derived the complete (compressible) NS starting from 
an action different from the previous works. Then it is found that the difference between 
f and f plays a crucial role and the velocity field in NS corresponds to f m . In below, we 
reproduce this derivation within the framework developed here. 



As is discussed in Ref. 



161 ]. we consider that the potential term, e (internal energy 



density) depends not only on p m but also on the other parameter s representing the internal 
thermodynamical state of the fluid element. In thermal equilibrium, s is identified as the 
specific entropy. Now, we write the variation as 



/-A 



J dt J d 3 Rj2p m0 E 



Ag ■ (DDx (R, t) + — VP ) + —6 x s 



1 



T 



711 



where T and P are defined in analogy to the usual thermodynamic relation as 

d{me j p^) 



+ 0(A 2 ), (77) 



T 



P = - 



Os 




d(l/p m ) \p 



(78) 
(79) 
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In the thermodynamic equilibrium, these quantities coincide with the usual temperature 
and pressure. 

Because we are considering dissipative processes, the entropy of a fluid element is not a 
conserved quantity. Then, to complete the variations, we need information about 5\s for 
the variation (1701) . In general, s is a scalar function of local quantities (does not depend 
explicitly on R) determined only by the properties of the fluid element. Therefore, the most 
general form of s depends only on the density p m and its derivatives, p m and Vp. m . Due to 
the vector nature of Vp m , this does contribute only to the higher order corrections of NS. 
Thus, the most general form of s for our purpose is expressed as s = s(p m ,p m ). 

Suppose that the hydrodynamic time scale Thyd = Pm/ Pm is much larger than that of 
the microscopic degrees of freedom inside a fluid element, r m j C . Note that the local thermal 
equilibrium should be perfectly achieved for the ideal fluid where s is constant (8\s = 0). 
Therefore, we expect that the specific entropy is expressed in powers of T mic /r hyd in such 
a way that the ideal fluid case is recovered in the vanishing limit of r m j C . We then write 
s = s eq + aiT mic p m / p m + a2{r mic f> m / p m ) 2 + ■ ■ • , where s eq is a constant for the variation and 
aj's are expansion coefficients. Thus the lowest order truncation gives 

6xs = 5x{g Tmic (Pm)Pm), (80) 

where g Tmic {Pm) is an arbitrary function of p m . For the stochastic variation, p m is interpreted 
as (D + D)p m /2. Note that 5\s is the virtual change of s associated with the variations and 
does not necessarily satisfy the thermodynamic principles such as 5\s > 0. 
Then the condition of 5\I = leads to 
d 



Pr, 



f + V(P-AiV-f ro )-(V-77V)f=0, (81) 



where rj is the shear viscosity. Here we used that v = f]/p m and p = 
—p m g(Pm)/ m (dT/dp m )s. One can see that the contribution from 5\s effectively changes 
pressure by /iV ■ f m . 

Equation ( 1HTT) contains the two different velocity fields, f m and f , but the latter can be 
eliminated by using Eq. f l30j) . As a result, we obtain 



|r + f m -V 

at 



f m + V(P-CV-f m )- V-(r/e* 



-(V-r/V) ( ^-Vln Pm ) =0, (82) 
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where e™ = {e™} is a (3 x 3) irreducible symmetric tensor, defined by 

e?i = dj + d t fj - |(V • f m )^-, (83) 

and ( = \i + 2?7/3. Identifying ( as the bulk viscosity, Eq. f|82|) is nothing but the Navier- 
Stokes equation, except for the last term on the left hand side. 

The last term is not only of second order for the magnitude of fluctuations (noise) v = 
v/Pm, but also of third order for the spatial derivative V. The Navier-Stokes equation 
does not contain such a term since, by construction, only the first and second order spatial 
derivative terms are maintained. In the case of rarefied gases, this corresponds to the first 
order truncation in the derivative ( Chapman- Enskog) expansion of one-particle distribution 
functions. 

On the other hand, in our case, the viscosity coefficients are directly related to the size 
of noises. When the effect of fluctuations is large, we cannot neglect the last term of Eq. 
( 1821 . This term comes from the difference between f and f m . Therefore, in such a case, 
hydrodynamics is described by using the two velocities as E q. (IHTl). This feature resembles 
that of the generalized hydrodynamics proposed by Brenner 
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22]. 



B. Gross-Pitaevskii Equation 



The Gross-Pitaevskii equation is a non-linear Schrodinger equation and has been used 



The derivation of the Gross- 



to study some aspects of the Bose-Einstein condensates 
Pitaevskii equation (GP) also has been studied from various points of view [24|. In particular, 
it is shown that GP can be derived in terms of SVM, extending the method of Sec. IIIIDI to 
a system of a many-body Lagrangian llj. In this subsection, we show that the formalism 
of SVM for a continuum medium can readily be used to derive GP in a direct way. 

We consider a system composed of identical bosons with mass m, under the external 
trapping potential ^(r). Supposing the interaction among bosons is local, we write the 
potential energy density of the system as 



PI 



Til 



(84) 



where p m = p rn (r, t) is the mass density of the boson and Uq represents the strength of 
interaction among them. To obtain the equation for the wave function (order parameter), 
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the kinetic term of the action should be replaced by the average of Dx and Dx as was done 
in Sec. IIIID1 Then, in the Lagrangian coordinate system, we can define the action of SVM 

as 



dt j d 3 Rp m0 (R) 



1 Dx ■ Dx + Dx ■ Dx 1 Tr/ . 
V x 

2 2 m 



2m : 



;U p m (x,t) 



(85) 



0. 



(86) 



Following the similar procedure in the previous subsection, SVM leads to 

(dt + f m • V)/;, - 2^ 2 V i (p- 1/2 V 2 v /p^) + -V 4 (v + % J 

my m J 

which should be solved together with the Fokker-Plank equation for p m , i.e., Eq. (168]) . 
This equations is same as that discussed in Ref. 25[ and known as the quantum Bernoulli 
equation. 

To convert Eq. (186]) into the corresponding Schrodinger-type equation, we follow the 
same procedure discussed in Sec. IIIIDl introducing the wave function (see Eqs. fp£?]) and 
(HT|) ) with a modification of Eq. fHTl) by 



, / Pm if) 

^ = \ — e , 
V m 



(87) 



because p m is the mass density. We again assumed that the flow is irrotational. The equation 
corresponding to Eq.f HHl) becomes 



id t t/j 



-u\/ 2 + -^—(V + U \ 
Imv 



(88) 



and setting v = h/(2m), we arrive at GP, 



ihd t ip 



2m 



\7 2 + V + U W 



(89) 



As before, from the definition of ip, \ip\ 2 = p m /m gives the particle distribution function. 

Note that the classical limit of GP is easily obtained by taking the vanishing limit of h 
(or v) in Eq.(|86l). 

(dt + f m • V)/;, = --V l V - ^V> m , (90) 
m m z 

together with Eq. (|68|) . Of course, these equations are same as those obtained from the 
action of Eq. ( 185]) applying the classical variational method. 
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The above argument is easily extended to include the external gauge fields via the min- 
imum coupling as described in Sec. IIVI The required action in the Lagrangian coordinate 
system is 



dt I d 3 TLp 



1 e 1 1 e 
v 2 + _ v . A v( x ) _ UoPm (y. t) _ _ ( 

2 cm m 2m z m 



(91) 



and the equation corresponding to Eq.fjoTj) becomes 

(d t + f m • V)f m - 2i. 2 V(p m 1/2 V 2 v^) + -V ( V + e0 + -U oPm 

m V m 



—d t A + — f m x (V x A). 

cm cm 



(92) 



The corresponding equation for the wave function is given by 



ihd t ip 



2^(-^V-^A) 2 +(y + e0 + W| 2 ) 



(93) 



where if) is given by Eq. flHTj) with the phase function given by Eq. fl58|) . and the noise is set 
as v = h/ (2m). 



C. Generalized Diffusion Equation and Physical Meaning of Correction to 

Navier-Stokes Equation 

In Sec. IV A| we found that NS derived using SVM involves the higher order correction 
term, which was neglected for deriving NS. In this section, we discuss the physical meaning 
of this correction term. 

For the sake of simplicity, we consider the case where there is no potential energy density 
in the action (R)3"|) . Physically, this situation will correspond to the free diffusion process. 
Then, from Eq. ( IT6]) . together with the continuity equation, we have, 

d t pm + V • (p m f m ) = 0, (94a) 

Pm(d t + f m • V)/ m * - ^djivpmdifj + vp m dj f^} - dj {{vp m )dj (vdi In p m )) = 0. 

i 

(94b) 

The first equation is the continuity equation as before, and the second one is obtained by 
SVM, that is, Eq. (1761) . In the discussion of Sec IV Al the third term on the left hand side 
of Eq. (194bl) was neglected as the higher order term of NS, but now we will keep it in this 
subsection. 
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This equation describes the generalized diffusion process. Suppose that there is a clear 
separation of time scales between p m and f m in the Lagrangian coordinate, and we observe 
the macroscopic dynamics of p m where the material derivative of f m is negligibly small. 
Then, from the second equation, we obtain 

f m = -^Vlnp m . (95) 

This corresponds to Fick's law of the diffusion equation. In fact, by substituting it into the 
first equation, we obtain the usual diffusion equation (see also Appendix IA1. 

dtPm = ^V 2 p m - (96) 

Differently from the diffusion equation, Eq. (I94bl) involves a memory effect through the 
time derivative of f m . In this form, we see that Eq. ( I94b|) is interpreted as the generalized 
diffusion equation, and it reproduces the usual diffusion equation in the Markov limit where 
the material derivative of f m is neglected. 

In other words, to reproduce the diffusion equation, we cannot neglect the third term on 
the left hand side of the second equation, which corresponds to the higher order correction 
term for NS. In this sense, this higher order term may have a special meaning and be an 
important correction even in NS. 



VI. CONCLUDING REMARKS 



In this paper, we discussed the application of the stochastic variational method to par- 
ticle systems and continuum mediums. The application to the particle systems has been 
studied from various points of view. First, we shortly summarized these results showing the 
derivation of the diffusion-type equation and the Schrodinger equation with and without the 
minimal gauge coupling. 

We further generalized of the idea of SVM to the calculations of continuum mediums, 
and showed that the Navier-Stokes, Gross-Pitaevskii and generalized diffusion equations are 
derived. 

As is known, the usual diffusion equation is inconsistent with underlying microscopic 
dynamics and should be modified 271, |28|. In fact, it is shown that the coarse-grained 
equation obtained from microscopic Hamiltonian using the projection operator method is not 
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the usual diffusion equation but the Maxwell-Cattaneo-Vernotte-type generalized diffusion 
equation. This equation has several advantages compared to the usual diffusion equation. 
For example, this generalized equation is consistent with an exact relation obtained from 
microscopic dynamics, while the diffusion equation is not 28|, |29|. Here, we proposed a 
new generalized diffusion equation in the framework of SVM. It is interesting to ask the 
relation between this generalized equation from SVM and the Maxwell-Cattaneo-Vernotte- 
type equation. 

We found that SVM predicts the higher order correction term to NS. This term becomes 
important when there is strong inhomogeneity of the mass density. In other words, we 
cannot recognize this effect from the study of incompressible fluids. The importance of the 
correction term becomes more clear in the discussion of the generalized diffusion equation. 
In order to reproduce the usual diffusion equation in the Markov limit, the same higher 
order term in the generalized diffusion equation plays an important role. In this sense, this 
term may be an important correction even to NS. 

There are several attempts to derive macroscopic equations for high density matters by 
assuming a Langevin dynamics for underlying microscopic degrees of freedom 30j. There, 
the interaction of the Langevin dynamics is assumed in advance, while it is determined 
by applying the variational principle in our approach. It is thus interesting to discuss the 
meaning of the correction term which we found, comparing with the behavior of such a 
dynamical density functional theory. More detailed study of this correction term is a future 
task |38( | . 

The advantage of SVM is that dissipative equations can be derived systematically once 
the action of the corresponding reversible dynamics is given. For example, in the physics 
of soft matter, dissipative equations are usually derived by applying Onsager's variational 



principle 32|, |33|. SVM may serve another convenient approach to discuss the dynamics of 



such a complex fluid. 

So far, we have assumed that the noise of SVM is a Gaussian white noise. It might be 
possible to generalize this discussion to the case of multiplicative noises. For the present 
work, we assume that the property of the noise is unchanged both for the forward and 
backward SDEs. This seems to be reasonable assumption because the properties of the noise 
is completely independent of the dynamics of macroscopic variables (Brownian particles) in 
this case of the additive noise. This is, however, not trivial for the multiplicative noise. If this 
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assumption is still applicable, the extention of the present argument to the multiplicative 
noise will be straightforward. 

As was mentioned in the introduction, the essence of the classical variational principle 
becomes more fundamental in the application of the approach to field theories and relativistic 



systems. As far as we know, SVM has not yet applied to such cases 34J. Studies of 
relativistic heavy ion collisions revealed that it will be of the fundamental importance to 
discuss dissipative mechanisms in relativistic systems. Thus, an extension of the present 
approach to relativistic field theoretical systems is expected to furnish an important clue 
to formulate such processes from the fundamental concept of coarse-graining procedure to 
arrive at the relativistic stochastic differential equation. However, the Lorentz covariant 



SDE still stays an open question [35]. 



This work was financially supported by CNPq, FAPERJ, CAPES and PRONEX. 

Appendix A: another limit to diffusion equation 

In Sec. IV CI the generalized diffusion equation is derived and the equation is reduced to 
the usual diffusion equation in the Markov limit. Interestingly, this usual diffusion equation 
is still reproduced when a special initial condition is chosen. 

Suppose that, as the initial condition of f m , we have 

v m = — i/Vlnp. (Al) 



Then we can easily check that both of the first and second equations of Eq. ( I94bl) are 
reduced to the same equation, 

d t Pm = vV 2 p m . (A2) 

That is, the generalized diffusion equation is reduced to the usual diffusion equation even 
when we use the special initial condition. The difference from Eq. (1961) is only the magnitude 
of the diffusion coefficient. 
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